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The masses, pseudoscalar and vector weak decay constants and electromagnetic 
form factors of light 5-wave mesons are studied in the framework of the relativis- 
tic quark model based on the quasipotential approach. We use the same model 
assumptions and parameters as in our previous investigations of heavy meson and 
baryon properties. The masses and wave functions of the ground state and radially 
excited tt, p, K, K* and (j) mesons, obtained by solving numerically the relativis- 
tic Schrodinger-like equation with the complete relativistic qq potential including 
both spin-independent and spin-dependent terms, are presented. Novel relativistic 
expressions for the weak decay constants of the pseudoscalar and vector mesons are 
derived. It is shown that the intermediate negative-energy quark states give signif- 
icant contributions which essentially decrease the decay constants bringing them in 
agreement with experimental data. The electromagnetic form factors of the pion, 
charged and neutral kaon are calculated in a broad range of the space-like momen- 
tum transfer. The corresponding charge radii are determined. All results agree well 
with available experimental data. 
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I. INTRODUCTION 



The theoretical investigation of the properties of light mesons such as tt, p, K, K* and 
is a longstanding problem which plays an important role in understanding the low-energy 
QCD. The description of these mesons within the constituent quark model presents addi- 
tional difficulties compared to heavy-light mesons and heavy quarkonia. In fact, due to the 
highly relativistic dynamics of light quarks, the v/c and l/m^ expansions are completely 
inapplicable in the case of light mesons and the QCD coupling constant as at the related 
scale /i is rather large. Moreover, the behaviour of as{^J^) in the infrared region is unknown 
and thus model dependent (exhibiting, e.g., freezing behaviour, etc.). The pseudoscalar 
mesons vr and K produce a special problem, since their small masses originate from their 
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Goldstone nature caused by the broken chiral symmetry. Therefore the rehable description 
of hght mesons requires the completely relativistic approach. It is well known that in the 
relativistic studies an important role is played by the Lorentz properties of the confining 
quark- ant iquark interaction. The comparison of theoretical predictions with experimental 
data can provide a valuable information on the form of the confining potential. Such in- 
formation is of great practical interest, because at present it is not possible to obtain the 
relativistic qq potential in the whole range of distances from the basic principles of QCD.^ 
Most of the main characteristics of light mesons are formed in the infrared (nonperturbative) 
region, thus providing an important insight in the low-energy properties of qq interaction. 
Thus investigation of both static (e.g., masses and decay constants) and dynamic (e.g., 
electroweak decay form factors) properties is of a significant importance. 

Many different theoretical approaches have been used for studying light mesons, which 
are based on the relativized quark model 0, the Dyson-Schwinger and Bethe-Salpeter equa- 
tions 0, 0], chiral quark models with spontaneous symmetry breaking (e.g. the Nambu- 
Jona-Lasinio model) (1| , the relativistic Hamilton dynamics p, Q , the finite-energy and 
light-cone [i| sum rules and lattice QCD [lOj]. Here we consider the possibility of investigat- 
ing light mesons on the basis of the three-dimensional relativistic wave equation with the 
QCD motivated potential. Our relativistic quark model was originally constructed for the 
investigation of hadrons with heavy quarks. It was successfully apphed for the calculation of 
their masses and various electroweak decays 11, ^ 3, JJ, 3|- ■'■^ these studies the heavy 



quark expansion has been used to simplify calculations. We determined all parameters of 
our model from few experimental observables (some masses and decay rates) and keep them 
fixed in all our subsequent calculations, thus ensuring its universality. While describing the 
properties of heavy-light mesons 2| , we treated the light quarks in a completely relativistic 



way. Recently this approach was applied for calculating the masses of light mesons [IC 
and light diquarks inside the heavy baryons Due to the phenomenological character of 
our model we cannot reveal the origin of the chiral symmetry breaking and thus the model 
cannot describe the chiral limit and the Goldstone nature of the pion. We consider the pion 
as the purely bound state of the quark and antiquark with fixed constituent masses. 

In this paper we extend our previous studies of light mesons and describe their elec- 
troweak properties such as the weak decay constants and electromagnetic form factors. The 
investigation of decay constants and form factors is an important issue since it provides a 
very sensitive test of the light meson wave functions and, thus, of the quark dynamics in 
a meson. It requires the completely relativistic consideration of the corresponding decay 
processes including account for the relativistic transformation of the meson wave functions. 
The comparison with the available large set of experimental data tests the model predictions 
in a broad momentum range and helps to discriminate between different model assumptions. 

The paper is organized as follows. In Sec. |H] we briefly describe our relativistic quark 
model, formulate our main assumptions and give the values of parameters. Then in Sec. IIIll 
we present our results for the light meson masses 16>|, for selfconsistency. There the proce- 
dure of constructing the completely relativistic local potential of the light quark interaction 
in a meson is described. The obtained potential is applied for calculating the light S'-wave 
meson masses and wave functions. In Sec. llVl the novel relativistic expressions for the weak 
decay constants of pseudoscalar and vector meson are derived. Special attention is paid 



^ Recent calculations of the nonperturbative qq potential in continuum Yang-Mills theory in Coulomb gauge 
can be fomid in Ref. 
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to including all possible intermediate quark states. It is argued that the negative-energy 
contributions play an essential role. The calculated decay constants are compared with 
other predictions and experimental data. The electromagnetic form factors of pseudoscalar 
mesons are studied in Sec. The relativistic expressions for these form factors are obtained 
which take into account the contributions of negative-energy quark states and relativistic 
transformations of the meson wave functions from the rest frame to the moving one. The 
calculated form factors are plotted in comparison with experimental data. The charged radii 
of the pion, charged and neutral kaon are also determined. Our conclusions are given in 
Sec. EH 



II. RELATIVISTIC QUARK MODEL 

In the quasipotential approach a meson is described by the wave function of the 
bound quark-antiquark state, which satisfies the quasipotential equation of the Schrodinger 
type till] 

where the relativistic reduced mass is 

_ ErE2 _ - {mj - mlf 

and E'l, E2 are given by 

= 2M ' = 2M ■ 

Here M = Ei + E2 is the meson mass, mi^2 are the quark masses, and p is their relative 
momentum. In the center-of-mass system the relative momentum squared on mass shell 
reads 

2 [JVP - (mi + m2)^][M^ - (mi - m^)^] 

The kernel V{p, q; M) in Eq. is the quasipotential operator of the quark-antiquark in- 
teraction. It is constructed with the help of the off-mass-shell scattering amplitude, projected 
onto the positive energy states. Constructing the quasipotential of the quark-antiquark in- 
teraction, we have assumed that the effective interaction is the sum of the usual one-gluon 
exchange term with the mixture of long-range vector and scalar linear confining potentials, 
where the vector confining potential contains the Pauli interaction. The quasipotential is 
then defined by ^ 

V{p, q; M) = Ui{p)u2{-pMp, q; M)ui{q)u2{-q), (5) 

with 

v(p, q; M) ^ v(p - q) = ^a,D,M7^i2 + KL(k)r^r2;^ + v;L(k), 



^ In our notation, where strong annihilation processes are neglected, antiparticles are described by usual 
spinors taking into account the proper quark charges. 
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where is the QCD couphng constant, Dni, is the gluon propagator in the Coulomb gauge 



An 



An 



D^' = D 



iO 



(6) 



and k = p — q; 7^ and u{p) are the Dirac matrices and spinors 



e{p) + m 



\ Hp) 



1 

crp 



e{p) + m 

with e{p) = ^p^ + m^. The effective long-range vector vertex is given by 



(7) 



%K 



r^(k) = 7^ + l^^v^vk"", 



(8) 



where k is the Pauli interaction constant characterizing the anomalous chromomagnetic 
moment of quarks. Vector and scalar confining potentials in the nonrelativistic limit reduce 
to 



= (1- 



reproducing 



Konf(r) = V^Ar) + K 



'conf ( 



Ar + S, 



(9) 



(10) 



illil. The 



where e is the mixing coefficient. 

All the model parameters have the same values as in our previous papers 
light constituent quark masses m„ = = 0.33 GeV, = 0.5 GeV and the parameters 
of the linear potential A = 0.18 GeV^ and B = —0.3 GeV have the usual values of quark 
models. The value of the mixing coefficient of vector and scalar confining potentials e = —1 
has been determined from the consideration of charmonium radiative decays Tl^. Finally, 
the universal Pauli interaction constant k = —1 has been fixed from the analysis of the fine 



splitting of heavy quarkonia Pj- states [ll|. In the literature it is widely discussed the 



't Hooft-like interaction between quarks induced by instantons [17[. This interaction can be 



partly described by introducing the quark anomalous chromomagnetic moment having an 
approximate value k = —0.744 (Diakonov jjlTi]). This value is of the same sign and order of 
magnitude as the Pauli constant k = — 1 in our model. Thus the Pauli term incorporates at 
least part of the instanton contribution to the qq interaction.^ 



III. LIGHT MESON MASSES 

The quasipotential (jSj) can be used for arbitrary quark masses. The substitution of the 
Dirac spinors (|7j) into (0) results in an extremely nonlocal potential in the configuration 



^ As is well-known, the instanton- induced 't Hooft interaction term breaks the axial [/yi(l)-syninietry, the 
violation of which is needed for describing the 1] ~ rj' mass splitting. We do not consider this issue here. 
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space. Clearly, it is very hard to deal with such potentials without any additional transfor- 
mations. In oder to simplify the relativistic qq potential, we make the following replacement 
in the Dirac spinors: 



ei,2(p) = Vmfo + p2^Ei,2 (11) 



(see the discussion of this point in [ij, |16|). This substitution makes the Fourier trans- 
formation of the potential (0) local. We also limit our consideration only to the S-wave 
states, which further simplifies our analysis, since all terms proportional to vanish as 
well as the spin-orbit ones. Thus we neglect the mixing of states with different values of 
L. Calculating the potential, we keep only operators quadratic in the relative momentum 
acting on Vcoub Korrf ^^"^ replace — > El 2 — "^1,2 i'^ higher order operators in accord with 
Eq. (jllll preserving the symmetry under the (1 2) exchange. 

The substitution works well for the confining part of the potential. However, it 
leads to a fictitious singularity 5^(r) at the origin arising from the one-gluon exchange 
part ( A Vcoui (?")), which is absent in the initial potential. Note that this singularity is not 
important if it is treated perturbatively. Since we are not using the expansion mv/c and are 
solving the quasipotential equation with the complete relativistic potential, an additional 
analysis is required. Such singular contributions emerge from the following terms 



[ei(?)(ei('?) + mi)ei{p){ei{p) + m^)] V2 

V^coui(k^), (12) 



[ei(g)ei(p)e2(g)e2(p)]^/ 



if we simply apply the replacement (fTT|) . However, the Fourier transforms of expressions ()12|) 
are less singular at r — >■ 0. To avoid such fictitious singularities we note that if the binding 
effects are taken into account, it is necessary to replace ei,2 -^1,2 — ''?i,2^, where V is 
the quark interaction potential and 771^2 = ""^2, 1/(^^1 + ^2)- At small distances r — > 0, 
the Coulomb singularity in V dominates and affords the correct asymptotic behaviour. 
Therefore, we replace ei^2 -^1,2 — ''?i,2Vcoui in the Fourier transforms of terms (fT^ (cf. 
[isl). We used the similar regularization of singularities in the analysis of heavy-light meson 
spectra 3|. Finally, we ignore the annihilation terms in the quark potential since they 



contribute only in the isoscalar channels and are suppressed in the ss vector channel 0. 
The resulting qq potential then reads 

V{r) = Vsi{r) + VsT,{r), (13) 

where the spin-independent potential for S-states (L^ = 0) has the form 



+ + + 

mim2 \ 

Ei + iTLi _^ -^1 + "^2 

El E2 

if i 

4 \Ei{Ei+mi 



{El + mi){E2 + m2 



+ T ( TTTT^ rAl^Coul 



1 



mi{Ei + mi) m2{E2 + m2) 




+ ^T^. r-(l + «:)(TT^ + TT^) AKL( 
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+ 



8mim2{Ei + mi)(i?2 + "^2 
and the spin-dependent potential is given by 

2 



AKo„f(^) 



(14) 



VsBir) 



3 El E2 



AVboui(r) + 1 — (1 + K 



2mi 



2mi 



Eo 



1712 



E2 + 777,2 



2mo 



S1S2 



(15) 



with 



V^li(r) 



_4a, 
3 r ' 
Vcoui(r) 



4a, 1 



3Ei + mir 



Vcoui(r) = Vcoui(r) 



4 a, 1 
3 Eir 



4 a, 1 
3-c;2 r 



Vl,2 



1,2) 



7772,1 



7r7i + 7r72 



(16) 



Here we put a, = aJu^) with /ii2 = 2mim2/ {mi + 7^2). We adopt for as(/x^) the simplest 
model with freezing |l9|, namely 

as(/i^) = - 



/3oln 



A2 



2 ' 



/3o 



11 - 3%, 



(17) 



where the background mass is Mb = 2.2A\J A = 0.95 GeV [19|], and A = 413 MeV was fixed 
from fitting the p mass. ^ We put the number of flavours nj = 2 for vr, p, K* and nj = 3 
for (j). As a result we obtain Cis{pli_d) = 0.730, ^^(pD = 0.711 and as(/i^g) = 0.731. 

The quasipotential equation ^ is solved numerically for the complete relativistic poten- 
tial (jl3j) which depends on the meson mass in a complicated highly nonlinear way. The 
obtained meson masses are presented in Table E] in comparison with experimental data [lOl 
and other theoretical results 0, S ■ This comparison exhibits a reasonably good overall 
agreement of our predictions with experimental mass values. Our results are also consistent 
with mass formulas derived using the finite-energy sum rules in QCD and with predictions 
of lattice QCD 10]. We consider such agreement to be quite successful, since in evaluating 
the meson masses we had at our disposal only one adjustable parameter A, which was fixed 
from fittin g th e p meson mass. All other parameters are kept the same as in our previous 
papers jll .11211. The obtained wave functions of the light mesons are used for the calculation 
of their decay constants and electromagnetic form factors in the following sections. 



IV. DECAY CONSTANTS 



The decay constants fp and fy of the pseudoscalar (P) and vector {V) mesons parame- 
terize the matrix elements of the weak current = qiJ^q2 = ^i7/i(l ~ 75)92 between the 



^ The definition (|17|l of can be smootlily matclied witli tlie ag used for heavy quarkonia at the scale 
about nic- 
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TABLE I: Masses of light 5'-wave mesons (in MeV) 



Meson 


btate 

25+1 r 




1 neory 






Experiment 
FDG I20J 


this work 






14] 


TT 


1 1 
VSq 


154 


150 


138 


140 


139.57 


P 


1 Dl 


776' 


770 


742 


785 


775.8(5) 


1 

TT 


^ So 


1292 


1300 




1331 


-| /~n~\ / -1 r\r\\ 

1300(100) 


/ 

P 




1 A Of 

140D 


1450 




1420 


1465(25) 


ft 

w 


3 Sq 


1788 


1880 




1826 


1812(14) 






1921 


2000 




1472 




K 


1 1 a 


482 


470 


497 


506 


493.677(16) 


K* 




897 


900 


936 


890 


891.66(26) 


K' 




1538 


1450 




1470 




K*' 


23^1 


1675 


1580 




1550 


1717(27) 


K" 


3^50 


2065 


2020 




1965 




K*" 




2156 


2110 




1588 




4> 




1038 


1020 


1072 


990 


1019.46(2) 


4>' 


23^1 


1698 


1690 




1472 


1680(20) 



fitted value 

corresponding meson and the vacuum. They are defined by 

(0|gi7^75g2|P(K)) = ifpK^, (18) 
(0|gi7'^g2|V(K,£)) = jyMye^ (19) 

where K is the meson momentum, and My are the polarization vector and mass of the 
vector meson. This matrix element can be expressed through the two-particle Bethe-Salpeter 
wave function in the quark loop integral (see Fig. H} 

(0|jf |M(K)) = / ^Tr{7,(l-75)^(M,p)}, (20) 

where the trace is taken over spin indices. Integration over p° in Eq. (j2U|) allows one to pass 
to the single-time wave function in the meson rest frame 

v^(M,p)= /^vl/(M,p). (21) 

This wave function contains both positive- and negative-energy quark states. Since in the 
quasipotential approach we use the single-time wave function \&a4'k(p) projected onto the 
positive-energy states it is necessary to include additional terms which account for the con- 
tributions of negative-energy intermediate states. The weak annihilation amplitude ()20|) is 
schematically presented in the left hand side of Fig. ^ The first diagram on the right hand 
side corresponds to the simple replacing of the single-time wave function ()21|) \E'(M, p) by the 
quasipotential one ^E'a/kIp)-^ The second and third diagrams account for negative-energy 



^ The contributions with the exchange by the effective interaction potential V which contain only positive- 
energy intermediate states are automatically accounted for by the wave function itself. 
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^(M,p) 





(1) 



(2) 



(3) 



(4) 



FIG. 1: Weak annihilation diagram of the light meson. Solid and bold lines denote the positive- and 
negative-energy part of the quark propagator, respectively. Dashed lines represent the interaction 
operator V. 



contributions to the first and second quark propagators, respectively. The last diagram 
corresponds to negative-energy contributions from both quark propagators. 
Thus in the quasipotential approach this decay amplitude has the form 



o|j^|m(k; 



^ Ml(pi) J'7w2(P2)^Mk(p) + 



\y (2vr)3 



Mi(pi)ri 



(27r)6 

V2U2{V2)y{V - p')^Af k(p) + (1^2 



+ 



M + ei(pO -e2(p'; 



(2vr) 



-ui(pi)ri 



M + ei(y) +e2(p') 



(22) 



where p^2 = K/2 ± p*^'); matrices Fi 2 denote the Dirac structure of the interaction poten- 
tial (0) for the first and second quark, respectively, and thus FiF2V^(p — p') = V(p — p'). 
The factor v2M follows from the normalization of the quasipotential wave function. The 
positive- and negative-energy projectors have standard definition 



aW(p) 



e(p) ± (m7° + 7°(7p)) 



The quasipotential wave function in the rest frame of the decaying meson \['a/(p) = ^a/o(p) 
can be expressed through a product of radial <l>Af (p), spin Xss' and colour 4>q\q2 wave functions 



<?1(?2- 



Now the decay constants can be presented in the following form 



JP,V — Jpy + JP,V + JP, 



(2+3) 



.(4) 



(23) 



(24) 



where the terms on the right hand side originate from the corresponding diagrams in Fig. 
and parameterize respective terms in Eq. (j22I)- In the literature 0, 21. 12^ 23l | usually only 
the first term is taken into account, since it provides the nonrelativistic limit, while other 
terms give only relativistic corrections and thus vanish in this limit. Such approximation 
can be justified for mesons containing heavy quarks. However, as it will be shown below, for 
light mesons other terms become equally important and their account is crucial for getting 
the results in agreement with experimental data. 

The matrix element (j^^ and thus the decay constants can be calculated in an arbitrary 
frame and from any component of the weak current. Such calculation can be most easily 
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performed in the rest frame of the decaying meson from the zero component of the current. 
The same results will be obtained from the vector component, however, this calculation is 
more cumbersome since here the rest frame cannot be used and, thus, it is important to 
take into account the relativistic transformation of the meson wave function from the rest 
frame to the moving one with the momentum K (see Eq. ()H4|1 below). It is also possible to 
perform calculations in the explicitly covariant way using methods proposed in '24'. 
The resulting expressions for decay constants are given by 



J py 




2ei(p) J 



ei{p) + mi][€2{p) +1712 




(25) 



.(2+3) 

J py 




ei(p)[ei(p) + ^1 



' e2{p) +1712 
\ 262 (p) , 

ti{p) + mi I ^ 
62 (p) + m2 J 



1/2 r 



M-6i(p) -62(p) 
_M+6i(p) -62(p) 



2) 



*p,y(p), 



(26) 



f(4) 

jpy 




'ei{p) +m{ 
. 261 (p) , 



, /62(p) +m2V'^ M- 61 (p) - 62 (p) 



2e2(p) 



M + 6i(p) +62(p) 



6i(p) + mi][e2{p) + m2] 



+ 



ei(p)e2(p) [ei(p) + rni][e2{p) + m2]_ 



*p,y(p), 



(27) 



with Ap = — 1 and Xy = 1/3. Here e is the mixing coefficient of scalar and vector confining 
potentials Q and the long-range anomalous chromomagnetic quark moment /t (jHl) is put 
equal to —1. Note that fp^^^ vanishes for pseudoscalar mesons with equal quark masses, 
such as the pion. The positive-energy contribution reproduces the previously known 



expressions for the decay constants j^|21|. The negative-energy contributions and (jlTj) 
are new and play a significant role for light mesons (see below). 

In the nonrelativistic limit p^/m^ — ^ the expression for decay constants gives the 
well-known formula 



NR 

py 



12 



M, 



1^ 



py 



py 



(28) 



where \&py(0) is the meson wave function at the origin r = 0. All other contributions vanish 
in the nonrelativistic limit. 

In Table 101 we present our predictions for the light meson decay constants calculated 
using the meson wave functions which were obtained as the numerical solutions of the 
quasipotential equation in Sec. 11111 The nonrelativistic values as well as the values 

of different contributions in Fig. [T] f^/'^'^'^^ fl ^ - (j77j) and the full relativistic results (EH) 



are given. In Table UTH we compare our results for the decay constants /m with predictions of 
other approaches ^ 2e A, 2s recent values from two- ^lO] and three-fiavour lattice QCD 
jiol l and available experimental data . It is clearly seen that the nonrelativistic predictions 
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TABLE II: Different contributions to the pseudoscalar and vector decay constants of light mesons 



(in MeV). The notations are taken according to Eqs. 1)24(1 and H28|). 

Constant fff fS^'^ + fj^' (ij^^ + 

1290 515 -391 -76% 124 

fx 783 353 -198 -56% 155 

fp 490 402 -183 -46% 219 

fx* 508 410 -174 -42% 236 

U 511 415 -170 -41% 245 



TABLE III: Pseudoscalar and vector decay constants of light mesons (in MeV). 
Constant this work [2j [3, J5] \7j Lattice \10} Lattice [26] Experiment [20] 





124 


180 


131 


219 


138 


126.6 ±6.4 


129.5 ±3.6 


130.7 ±0.1 ±0.36 


fK 


155 


232 


155 


238 


160 


152.0 ±6.1 


156.6 ±3.7 


159.8 ± 1.4 ±0.44 


fp 


219 


220 


207 




238 


239.4 ± 7.3 




220 ± 2* 


fx* 


236 


267 


241 




241 


255.5 ±6.5 




230 ± 8t 


fs 


245 


336 


259 






270.8 ±6.5 




229 ± 3^ 



* derived from the experimental value for TpO^e+e-- 

^ derived from the experimental value for the ratio TT—tK*iyr f^r^pi/r ^^d the fp value, 
derived from the experimental value for F^^g+e- • 



are significantly overestimating all decay constants, especially for the pion (almost by a 
factor of 10). The account of the part of relativistic corrections by keeping in Eq. ()24)1 
only the first term /|}^ fl25)l . which is usually used for semirelativistic calculations, does not 
dramatically improve the situation. The disagreement is still large. This is connected with 
the anomalously small masses of light pseudoscalar mesons exhibiting their chiral nature. In 



the semirelativistic quark model |21| the pseudoscalar meson mass is replaced by the so- 



called mock mass Mp, which is equal to the mean total energy of free_ quarks in a meson, and 
with our wave functions: = 2{eq{p)) ^ 1070 MeV (~ 8M^) and Mk = (eg(p)) ± (e,(p)) 
1232 MeV (~ 2.5Mk)- Such replacement gives fp^ values which are still ~ 1.4 times larger 
than experimental ones (cf. 0). As we see from Table ITTl in the quasipotential approach it 
is not justified to neglect contributions of the negative energy intermediate states for light 
meson decay constants. Indeed, the values of Jm^^^ + /m"* large and negative (reaching 
—76% of fl^^ for the pion) thus compensating the overestimation of decay constants by 
the positive-energy contribution fl^\ This is the consequence of the smallness of the light 
pseudoscalar meson masses compared to the energies of their constituents. The negative- 
energy contributions ()26j) . ()27|) are proportional to the ratio of the meson binding energy 
M — ei(p) — e2{p) to its mass. For mesons with heavy quarks this factor leads to the 
suppression of negative-energy contributions since the binding energies are small on the 
heavy meson mass scale. This results in the dominance of the positive-energy term since 
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FIG. 2: Lowest order vertex function T^^^ corresponding to Eq. I\32\i. Photon interaction with one 
quark is shown. 



the negative-energy terms give only l/mq contributions {itlq is the heavy quark mass).^ On 
the other hand, for hght mesons, especially for the pion and kaon, the binding energies are 
large on the light meson mass scale and, thus, such factor gives no suppression. Taking the 
complete relativistic expression for decay constants /m (EH) brings theoretical predictions 
in good agreement with available experimental data. 

The comparison of our values of the decay constants with other predictions in Table lllll 
indicate that they are competitive even with the results of more sophisticated approaches 
j^, which are based on the Dyson-Schwinger and Bethe-Salpeter equations. On the other 
hand our model is more selfconsistent than some other approaches 0, S 2^1- We 



calculate the meson wave functions by solving the quasipotential equation in contrast to 
the models based on the relativistic Hamilton dynamics p, Q where various ad hoc wave 
function parameterizations are employed. We also do not need to introduce the mock meson 
mass 0, Uli 122] and to substitute it for the light meson mass as it was discussed above. 



V. ELECTROMAGNETIC FORM FACTORS 

The elastic matrix element of the electromagnetic current between the initial and final 
pseudoscalar meson states is parameterized by the form factor Fp{Q'^) 

{M{Pf)\J^\M{Pj)) = Fp{Q'){Pj + Pf)^, (29) 

where = -{Pp-Pif- 

In the quasipotential approach such matrix element has the form [l^] 

{M{Pf)\J,\M{Pi)) = I ^^^^r,,p^(p)r^(p,q)vl/^,p^(q), (30) 

where r^(p, q) is the two-particle vertex function and are the meson wave functions 
projected onto the positive energy states of quarks and boosted to the moving reference 
frame. The contributions to T come from Figs. |21 and El The term T^"^^ includes contributions 



^ For the heavy- heavy Be meson (cb) these negative-energy corrections will be of order /c^ and thus very 
small. The influence of the negative-energy contributions f'^'^^'^'^ on the decay constants of heavy-light 
B and D mesons will be considered elsewhere. 
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92 



FIG. 3: Vertex function T^'^^ corresponding to Eq. Dashed lines represent the interaction 

operator V. Bold lines denote the negative-energy part of the quark propagator. As on Fig. [21 
photon interaction with one quark is shown. 

from the negative-energy quark states. Note that the form of the relativistic corrections 
resulting from the vertex function T^"^^ explicitly depends on the Lorentz structure of the 
gg-interaction. Thus the vertex function is given by 



where 
and 



r^(p,q) = r«(p,q) + rf (p,q) + ---, 
rW(P,q) = eiMi(pi)7^Mi(gi)(27r)3(5(p2 - qa) + (1 2), 



r!f)(p,q) = ei^i(pi)^2(p2)^V(p2-q2)-7^U^^7?7i. 



+7lM 



, 7i°V(P2 - q2) j«i(gi)«2(g2) + (1 2) 



eiik 



Here ei,2 are the quark charges, ki = pi — A; k'^^ = qi + A; A = Pp — P/; 



(_) ^ e{p) - (m7° + 7°(7p)) 

Hp) 



e{p) = ypM-m^ 



and 



Pi,2 = ±j:^^hpf)p' 

i=i 

91,2 = el,2(9)§±Er^«(P/)9^ 

2 = 1 



where n^*-* are three four-vectors given by 



E= Vp' + M2, 2,j = 1,2,3, 



(31) 
(32) 



(33) 



Pi = {Ej, P/) and Pp = {Ep, Pp) are four-momenta of the initial and final mesons. 
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It is important to note that the wave functions entering the current matrix element ()3U|) 
cannot be both in the rest frame. In the initial meson rest frame, the final meson is moving 
with the recoil momentum A. The wave function of the moving meson \I/ma is connected 
with the wave function in the rest frame \1/mo = ^a/ by the transformation j27|| 



vI^ma(p) = D\'\R^jDr{R'i:j^Moi:p) 



1/2, 



(34) 



where is the Wigner rotation, is the Lorentz boost from the rest frame to a moving 
one, and the rotation matrix D^/'^{R) in the spinor representation is given by 




S-\vi,2)S{A)S{p), 



(35) 



where 



S{p) 



i{p) + m 
2m 



1 + 



e{p) + m J 

is the usual Lorentz transformation matrix of the Dirac spinor. 

To calculate the matrix element ()29|1 of the electromagnetic current between the pseu- 
doscalar meson states we substitute the vertex functions T^^) ^ and T^^) (jHH)) in Eq. (jHUl) 
and take into account the wave function transformation (j34p . Then we use the 6 function in 
r'-^^ to perform one of the integrations in the matrix element ()3()|1 . For the contribution of 
r*^^) we use instead the quasipotential equation to replace the integral of the product of the 
interaction potential and the bound state wave function by the product of the correspond- 
ing binding energy and the wave function. To simplify the calculation we explicitly use the 
value K = — 1 for the long-range anomalous chromomagnetic quark moment dHJ. However, 
as previously we keep the dependence on the mixing parameter e of the vector and scalar 
confining potentials Q- As a result we get the following expression for the electromagnetic 
form factor of the pseudoscalar meson: 



+ eFP'{Q') + {l-e)FP''{Q' 



(36) 



d^p - ( 262 (p) . 



(27r)3 
pA 



Xp) + mi 



E + M ;^ei(p+A)+mi 
ei{p + A) - ei{p) 



ei(p + A) +ei(p) 



2./ei(p + A)ei(p) 



Pt 



2^ei(p + A)ei(p)(ei(p)+mi) 2^/ei(p + A)ei(p) ^i{p) + 
( 1 1 



^Af(p) + (1 



(37) 



.(2)5 



2^EM ( f <fp 



2VJ^M f 



ei{p + A) - ei{p) + 2mi 



262 (p) 
E + M' 



ei {p) + mi ei (p + A) + mi 
^ei(p + A)+mi 2ei(p + A) 
pA 



. 2^ei(p + A)ei(p) 
ti{p + A) + mi 



Pr 



2jei(p + A)ei(p)(ei(p) +mi 
1 1 \' 



+ 



2w'ei(p + A)ei(p) ei(p) + \ei{p) + mi t2{p) 



m2 



14 



X 



ei(p + A) -ei(p) 



ei(p + A)[ei(p + A) + ei(p)] 



[M - ei(p) - e2(p)]^A/(p) + (1^2) 



(38) 



,(2)y 



X 



(27r)3 
ei(p) - mi 



Af P 



+ 



2e2(p) 



ei(p)+mi ei(p + A)+mi 



E + M y^ei(p + A)+mi 2ei(p + A) 
pA 



.2^ti{p + A)ei(p) 2^ei(p + A)ei(p)(ei(p) + 
ei(p + A) + mi p|, / 1 1 



+ 



X 



2^ei(p + A)ei(p) + "^i Vei(p) + "^i e2(p) 
ei(p + A) -ei(p) 



m2 



ei(p + A)[ei(p + A) + ei(p)] 



[M - ei(p) - e2(p)]*M(p) + (1 2) , (39) 



where p^^^^'^'^ are contributions from scalar (vector) confining potentials and pr = p^ — 
(pA)^/A^, E = ^ M"^ + A^. As previously, we put e = —\ for further numerical calcula- 
tions. It is important to note that the above expressions for the electromagnetic form factor 
of the positively-charged pseudoscalar meson exactly satisfy the normalization condition 



i^p(O) 



(40) 



following from the electric charge conservation. 

Now we can use the wave functions of the pseudoscalar light mesons (vr, K)^ found 
in Sec. IIIIl for the numerical calculation of their electromagnetic form factors Fp{Cf' 



m 



the space-like region Q > 0. The results of such calculations for the charged pion are 
shown in Figs. |^ {F^{Q'^)) and |S1 (Q^F^((5^)) in comparison with experimental data from 
Refs. (isl . liil 130^ . Good agreement with data both in low and high regions is found, 
including recent JLab data (3fl| which are plotted with crosses. It is clearly seen from Fig. El 
that the calculated pion form factor at high exhibits the asymptotic behaviour F^^IQ"^) ~ 
<^s{Q'^)/Q^ predicted by the quark counting rule [lH and perturbative QCD |3^. O ur results 
for the pion form factor can also be compared with QCD based calculations j33| and with 
recent parameterizations 



35| which arise from the constraints of analyticity and unitarity. 
The latter form factor models are based on the vector meson dominance and include a 
pattern of radial excitations expected from dual resonance models [3^. The consistency of 
our results with such parameterizations (cf. Fig. |3] with Fig. 2 of Ref. [s^) just means the 
manifestation of the quark-hadron duality. Finally, our predictions agree fairly well with 
recent lattice computations of the pion form factor |i3&.i37i]. The corresponding plots for the 
charged kaon form factor are given in Figs. IHl and [7| in comparison with experimental data 
from Refs. [i^ 3^|, which are available only for the low region. Again good agreement 
with experimental data is found. On Fig. [7|we also plot the neutral kaon form factor by the 
dashed hne. 

The mean-squared charge radius of the pseudoscalar meson (P = tt, K) is defined by 



-6 



dFpjQ' 
dQ2 



(41) 



The calculated values of the charge radii of light pseudoscalar mesons are given in Table |IV| 
in comparison with predictions of other approaches 0,0 EE 22 1 and experimental data |,2Qi] . 
An overall good agreement with experimental data is found. 
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FIG. 4: The charged pion form factor squared in comparison with experimental data from Refs. [21 
(open circles), |29| (solid squares) and |30|] (crosses). 



TABLE IV: Charge radii of light pseudoscalar mesons. 



charge radii 


this work 


\2] 


[25] 


m 


Lattice [37] 


Experiment [20] 


V{r''). (fm) 


0.66 


0.66 


0.67 


0.63 


0.63 ±0.1 


0.672±0.08 


V{r')Ki (fm) 


0.57 


0.59 


0.62 


0.60 




0.560±0.031 


{r')KO (fm^) 


-0.072 


-0.09 


-0.086 


-0.062 




-0.076±0.018 
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FIG. 5: times charged pion form factor in comparison with experimental data from Refs. ^ 
(open circles), [i^ (solid squares) and js^] (crosses). 
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g2 (GeV2) 

FIG. 6: The charged kaon form factor squared in comparison with experimental data from Refs. [si 
(open circles) and [s^ (solid squares). 



VI. CONCLUSIONS 



The relativistic quark model, which has been previously developed and successfully used 
for the comprehensive investigation of different properties of heavy and heavy-light hadrons, 
was applied here for calculating the masses, weak decay constants and electromagnetic form 
factors of the light mesons. The main assumptions and parameters of the model (such as the 
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FIG. 7: times the charged kaon (sohd hne) and neutral kaon (dashed hne) form factors. 



Lorentz structure and parameters of the confining potential and quark masses) were kept 
the same as in previous studies. The only change we made, is the necessary modification 
of the running coupling constant in the infrared region. Following Ref. we chose 

the simplest model with freezing (jl7|) . Therefore only one additional parameter A was intro- 
duced and it was fixed from fitting the p meson mass. We constructed the local relativistic 
quasipotential for the light quarks using the replacement (fTTj). which was previously tested 
on the heavy-light mesons. The resulting relativistic potential ()13p depends on the meson 
mass in a complicated nonlinear way. Solving numerically the quasipotential equation 
we got masses of the ground-state and radially-excited light mesons in a reasonably good 
overall agreement with experimental data. Even the masses of the pseudoscalar tt and K 
mesons are well reproduced. This is a nontrivial result, since we use the constituent quark 
masses in our description and thus the chiral symmetry is explicitly broken from the very 
beginning. We determined the light meson wave functions and used them for studying their 
electroweak properties. 

First the weak decay constants of pseudoscalar and vector mesons were investigated. 
It was argued that both positive- and negative-energy parts of the quark propagators in 
the weak annihilation loop should be taken into account. Usually in the semirelativistic 
quark model [2? |21|, l22|] only the positive-energy contributions are kept. This approximation 
requires to replace in the expression for the pseudoscalar decay constant (P^j) the meson 
mass by the so-called mock meson mass, which is considerably larger, in order not to get 
the significant overestimate of the decay constants. We showed that the negative-energy 
contributions to the light meson pseudoscalar decay constants are large and negative. Their 
account brings theoretical predictions (with the physical meson masses) in good agreement 
with available experimental data. 

Next we studied the electromagnetic form factor of the pseudoscalar mesons. The corre- 
sponding matrix element of the electromagnetic current was calculated using the quasipoten- 
tial approach. The additional contributions of the intermediate negative-energy states ()33|) 
were taken into account as well as the transformation of the meson wave function from the 
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rest frame to a moving one ()34|). As a result the relativistic expression for the electromag- 
netic form factor was obtained. We then calculated the pion, charged and neutral kaon form 
factors in the space-like region. Good agreement with available experimental data both in 
small and large regions were found. At large momentum transfer this form factor tends 
to reproduce the power- law behaviour predicted by perturtbative QCD The calculated 
charge radii of light pseudoscalar mesons are in good agreement with experiment. 

In conclusion, we found that the obtained results are quite competitive with the predic- 
tions of other approaches 0, SB IS 0110, 25, 26, 31,^] including more sophisticated ones, 
which were specially developed for treating light mesons. 
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